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Abstract- This paper discusses RCS (Radar Cross Section) 
predicting software code based on high-frequency technique, 
Physical Optics (PO) and Physical Theory of Diffraction (PTD) 
which are powerful asymptotic technique for the surface and 
edge currents integral in the computation of RCS for large and 
complex object. By distinguishing between illuminated and un- 
i 11 unlimited areas of target for multiple scattering fields, we can 
obtain total polarimetric scattering fields of an observation point. 
Target model is available for not only isolated one but also wide 
area terrain surface such as ground or sea, with various forms of 
vegetation, tree, or building. In order to keep high shadowing 
accuracy, an adaptive subdivision scheme algorithm, which is 
depending on the situation of incident-, reflection-, scattering- 
shadowing and multiple reflections, is applied to processing of 
this distinction. Since both of co- and cross-polarized RCS fields 
are calculated, the results are useful and effective for studying 
full polarimetric radar image processing or remote-sensing image 
decomposition analysis. In this paper we show description of 
scattering theory for the developed code to calculate full 
polarimetric RCS and several calculated RCS examples including 
Doppler spectrum simulation and measurement. 

Keywords- Polarimetric Radar Cross Section (RCS); Physical 
Optics (PO); Physical Theory of Diffraction (PTD); Adaptive 
Shadowing; Systematic Prediction Code; Doppler spectrum 

I. INTRODUCTION 

The significance of the RCS in the outcome of radar 
engagements makes its prediction an important problem in 
modern systems such as polarimetric Synthetic Aperture 
Radar (SAR) on space/air-borne for remote-sensing or 
collision avoidance radar for Intelligent Transfer System (ITS). 
Many useful techniques have already been developed to 
provide realistic results. They can be based on exact methods 
(IE: integral equations, MoM: method of moments, FDTD: 
finite-difference time-domain, FEM: finite element method 
and so on [1, 2]) or high-frequency approximate approaches. 
One can refer to commercial or non-commercial software 
package such as WIPL-D, FEKO, XFDTD, XPATCH [3] and 
so on. Generally, most numerical methods are exact but have 
commonly a drawback to be limited to electrically small size 
objects. 

In this paper, we discuss the systematic software code, 
which is developed by using typical high-frequency method, 
namely the Physical Optics (PO) and Physical Theory of 
Diffraction (PTD), in order to predict the RCS of complex and 
large targets on or not on wide area such as ground. The PTD 
is an integrative powerful technique in which the PO uniform 
current on a discontinuous surface of target is refined by the 
additional non-uniform component due to the presence of the 
discontinuity especially such as an edge. The developed code 
has minimum computer resource requirements and provides 
convenient run-times. 



It utilizes the scientific and engineering computational 
features of MATLAB and its Graphical User Interface (GUI) 
functions to provide an error-free encoding of input parameters 
and efficient calculation. The developed code, named ACES 
(Accurate Calculation of Electromagnetic Scattering), provides 
many powerful computational capabilities, including the 
shadowing, the multiple reflections, the approximate effects of 
the ground and the effects of materials and coatings. Output 
result of RCS is not only co-polarized but also cross-polarized 
field. Four powerful engines are developed and embedded in 
the software. They are Incident Shadowing Engine (ISE), 
Reflection Shadowing Engine (RSE), Scattering Shadowing 
Engine (SSE) and Multiple Reflection Engine (MRE). By 
dynamically re-segmenting the triangular facets according to 
the situation of shadowing and reflection at the specific angle 
in the computation process, we can obtain good agreements 
with measurement. 

By employing PTD, the effect of current along edge is also 
considered in the code. Induced surface current representation 
can be related to the dynamic re-segmentation of the triangular 
facets which is depending on the situation of incident, 
reflection, scattering shadowing and multiple reflections. A 
user-friendly GUI design of materials database is also 
developed. A user updateable database of materials and 
coatings can be applied to models in one or multiple layers, and 
the computation of their effects on the RCS of models. The 
user can easily apply the information of material to the 
computational model and get the information of material which 
is applied on the computational model interactively. 

II. FORMULATION PO AND PTD BASED ON ELECTROMAGNETIC 

CURRENTS 

The Physical Optics (PO) approximation is one of the most 
convenient well known RCS prediction methods for an 
arbitrary 3-dimensional target [4], The induced currents are 
integrated over the illuminated portions of the target surface to 
obtain the scattered far field, while setting the current to zero 
over the shadowed portions. Thus, the electromagnetic surface 
currents are given as J s (r) = n x H l ,M s (r) = E l x n for 
illuminated region and zero for shadow region, where n is the 
unit normal vector of the surface directed outward and// 1 ,E l 
are the incident magnetic and electric field vector respectively. 
The electric and magnetic surface currents, J s (r ) and M s (r ), 
are then used in the radiation integrals to compute the scattered 
far-field from the target. PO is a high-frequency approximation 
method that gives best results for electrically large bodies 
(L > 5/1,2: wavelength) and is most accurate in the specular 
directions. If the source illuminating the target is at a far 
enough distance, then the incident field can be regarded as a 
plane wave. PO surface current is integrated over the 
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illuminated surface of the scatterer. The scattered field 
component is found by calculating the following integral 
equation [5], 



-jk 



exp(-jkR) 



4nR 



■ {ZoJs -h,<p±M s - Vfl} expQ'fcr' ■ \ R )ds (1) 
J s 

where k = 2n/A, A is the wavelength of the incident wave, r' 
is position vector for the source/ s (r ) and M s (r ), i fi and R are 
unit vector and the distance to observation from origin, 
%o = Vl I o/ £ o i s f ree space impedance, (i e .i^ ) are the unit 
vectors of the spherical coordinate system, and (Eg ,E?°) are 
the scattered electric field components which are evaluated 
over a fit region S in the spherical coordinates. Equation (1) is 
for co-polarized field and by exchanging i e and i^ , cross- 
polarized field can be obtained. 

For a case of non-metal flat plate, we can more explicitly 
derive a following reduced expression from the proposing 
Formula (1). 

Eg P H = -Jk^^jr 1 D e4 / s exp{-;fc(fc< -fe s ) r } dS' (2) 

D g = ^{(l + fi||)cos0 - (1 - fi||)cos0 o }cos(0 o - 0) 
-E 2 {(l - R L ) cos 6 cos 9 - (1 + R ± )} sin(0 o - 0), 

D = EjQ. +i?|,)- (1 -fl||)cos0 o cos0}sin(0 o - <p) 
+E 2 {(l -RJcosBq + (l + i? ± )cos0}cos(0 o -0), 

where E x and E 2 are amplitudes of the incident wave given as 
E l = \E-t\g + E 2 i ( p}exp(jkk i ■ r), fe'and fc s are wave-number 
vector of incident and scattered wave, and i? B and R L are 
parallel and perpendicular coefficients of reflected wave, 
respectively, as discussed below. Judgment for lit or shadow 
region can be processed by using formula k l • n < , which is 
an illuminated case. The incident angle (0 o ,0 o )and the 
observation angle (6,(p) are independent each other. For 6- 
polarized (TM-wave) or 0-polarized (TE-wave), amplitude of 
incident wave E 2 or E 1 becomes to be zero, respectively. In 
case of E 1 = E 2 , polarization of scattered field becomes to be 
slanted 45 degree fro m the unit vector i e and ij, . 

When a body is illuminated by a plane wave, current is 
induced on its surface. Around the edges of the body, this 
current has more complicated behavior, due to the diffraction 
of the incident field. The edge diffracted field appears to come 
froma non-uniform line source located at the edge. Concept of 
PTD method is an idea to consider this non-uniform current 
into PO additionally. According to the PTD or the method of 
equivalent currents, any finite current distribution yields a 
finite result for the far-zone diffracted field. The diffracted field 
E d due to an edge discontinuity C is given for far zones by the 
radiation integral as eqivalent edge currents [6, 7], 

E d = -jk j (Z i r xi r xi t If- 0')+irXi £ Mf u (I')} 

exp {—jkr) 



Anr 



dV 



(3) 



where, as shown in Figure 1, i t is tangential unit vector of the 
edge, If (!') and Mf u (l')are equivalent electromagnetic edge 



currents along edge length 1' , and r is distance between 
observation and the edge of target. Then the total scattered 
fields of PTD can be represented as E s = E po + E d . 




Figure 1: Diffraction by a wedge and equivalent current along the edge. 

III. MULTIPLE REFLECTION SCHEME FOR CURVED IMPEDANCE 

SURFACE 

Multiple reflection is considered that a reflected plane wave, 
which illuminates to an impedance plane #m, again incidents to 
an impedance plane #m+l. Regarding as subscript (m) means 
plate number and R^(m) l R L (m) are reflection coefficient for 
parallel 6- (TM-wave) and perpendicular <p -polarization (TE- 
wave), respectively, of plate #m, relating expressions are 



E \m) ~ P lm) T (m) + Q(m) B (m)' 

Elm) = Si, (m)P{ m) n m ) + R L(.m)Ql m) B Gn) 

Km) Xn (m) Km) Xn (m) 



B im) ~ 



sin# 



(m) 



sin9 



Cm) 



T (m + 1) ~ B (m) X Km)' 



\m) 



Km) 2n (m)[ n (m) " fc ( m )] 



(4a) 
(4b) 

(4c) 

(4d) 
(4e) 



where subscripts t,rare meaning to incident and reflected 
wave, respectively, P i,r ,Q i,r are amplitude and phase variable 
of geometrical wave. Expressing wave-number unit vector of 
incident and reflected wave are k l , k r (or k L , k r in Figure 2 and 
3) respectively, then a vector B vertical to incident plane is 
defined as B = k l x n/sin# = k r xn/sinB , and 8 Q is an 
angle between k l and n as shown in Figure 2. The vector k l 
andfc r are determined by the low of Snell, k r = k l — 2n(n • 
k l ) . A vector T iJ can be defined as T Ur = B x k ir . 




Figure2: Geometrical wave reflection from curved impedance surface. 
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flat plate scatterei 



Figure 3: Multiple reflections among impedance surfaces. 

Thus, wave number vector of reflected wave at a plane #m, 
k r (m y equals to kl m+1 , which is one of incident wave at a plane 

#m+l. In general B (m) and B (m+1) , T l £ n) and T^ +1) are not 
equal, as shown in Figure 3. So, if relative geometrical position 
relation of a plane #m and #m+l is once obtained, explicit 
expression can be expanded, especially normal vector relation 
of each plane is most important. In this situation, we can get a 
following relation [8], 



B 
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IV. REFLECTION AND REFRACTION COEFFICIENTS OF MULTI- 
LAYERED DIELECTRIC SHEETS 

In this section, we discuss reflection and refraction from 
multiple layered dielectric flat sheets, as shown in Figure 4. For 
scattering by a target which has dielectric portions, introducing 
boundary condition of surface impedance makes easy to 
formulate electromagnetic relations. On the other hand, we can 
also obtain exact general formulation for this model such a 
multiple layered dielectric sheets which is infinitely flat. As a 
result, the coefficients of reflection and refraction obtained 
from the solution of wave equation are exactly formulated by 
characteristic matrix of dielectric permittivity, magnetic per- 
meability and thickness of each dielectric sheet. And we can 
calculate the coefficients from metal backed sheets by infinite 
permittivity of a sheet and dielectric loss by considering 
complex permittivity [9]. 
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Figure 4: Multiple reflections between impedance flat planes. 



Now as shown Figure 4, N-layered sheets with infinite 
length to y-axis are overlapped along z-axis. Incident angle is 
8 i=1 measured from normal z-axis of a first sheet, s t and fi t 
(i = 1 ~ N) are permittivity and permeability, respectively, 
which are function of variable z. In case TE-wave, electric field 
vector is perpendicular to incident plane, components y, z of 
the field are 0, and then only a group of (E X ,H ,H Z ) remains. 
Concentrating to electrical field component, E x , Maxwell's 

i -• . dE x . dE x 

equation reduces to 3-equations, \- jcj/xHy = , 

dH z dtly 

jo)fj.H z = 0, j(i)£E x = . Assuming permeability 

3 (fHy ) (jli 3 Hy 

depends on u = u(z) and noting = H v \- u — Idz , 

'" " dz dz 



dz 



we can obtain following basic equation 



2 d 2 ) d[i d(]ogn)dE x 

dy 2 dz 2 ) y dz dz dz 



(6) 



where a) 2 £p. = n 2 k 2 . After some calculation, for example 
using variable separation method, following characteristic 
matrix for a dielectric sheet is obtained [9]. 



M^,cp) 



cos<p — —sincp 
,j!;sin<p cos q> / 



(7a) 



{ J s r /fj. r cos 8 ■■ TE-wave 
. (7b) 
Jn r /£ r cos9 ■■ TM-wave 

Letting the thickness of each sheet is z t (i = 1 — N), each 
characteristic matrix is recurrently obtained asi = 1 : G = 
M x G(z x ),i = 2 ■■ G(z t ) = M 2 G(z 2 ), ••• . Therefore, the matrix 
of the sheets i = 1 to N is 

G =M r G{z l ) = M 1 M 2 G(z 2 ) = ■■■ 

= M 1 M 2 --M N G(z N ) = M{z N )G(z N ). (8) 

The relation between the matrix elements and coefficients 
can be derived as follow. Regarding incident angle of first 
(t = 1) and last boundary as 6 X and 8 N , and amplitude of 
incident, reflected and refracted wave as A i ,A r , and A t , 
reflection and refraction coefficient are defined to 7? = A r / 'A t 
and T = A t /A it respectively. As m^ are the 4 elements of 
matrix M, and considering the constraint condition A t + A r = 
(m n + m 12 ^ N )A t and ^(4; - A r ) = fjn 21 + m 22 % N )A t , the 
coefficients are lastly obtained as 

A r m 11 ^ 1 -m 12 ^ 1 ^ N -m 21 -m 22 ^ N 
R = T" = : r^ — : : r~ , (9a) 



*t = ^£i 

A l m xx S\ +m 12 ^^ N +m 21 +m 22 ^ N ' 



(9b) 



For example N = 3, let first and third sheet are free space 
and permittivity and permeability of second (t = 2) sheet are 
£ 2 ,\i 2 , respectively. For this case, the characteristic matrix is 
calculated as 



m n = m 22 = cos(p,m 12 =;' — sin(p,m 21 = ;'f 2 sin cp, 

S2 



q> = n 2 kz 2 cos t9 2 , !; 2 = n 2 cos9 2 . 



(10) 



Therefore putting i = 3 into previous equation, the coefficients 
are obtained as 
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T = 
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R u - 
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£2 ~?3 "2 cos6, 2 -n 3 cos0 3 
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2fi 



n 2 cos# 2 + n 3 cos# 3 ' 
2n 1 cos8 1 
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^23 — 



fi + £2 ^icosffj + n 2 cos# 2 

2^2 



2n 2 cos8 2 
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(Ha) 
(lib) 
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(lid) 
(lie) 
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where n 12 = 7 £ r(i),r(2)/M r (i),r(2) forTE-wave. As expected, 
these coefficients are a single dielectric sheet and perfectly 
agree with Fresnel formula. 
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(b)Wiih a like taiget: N=7 

Figure 5: Multiple reflection application for wall-through-radar. 




Figure 6: Reflection (Left) and transmission (Right) coefficients with no 
high density dielectric (metal- like target) sheet. 




Figure 7: Reflection (L) and transmission (R) coefficients with high density 
dielectric (metal-liketarget) sheet on the same condition of Figure 6. 



Next we show an interesting and simple example for a like 
wall-through-radar indicated in Figure 5 configuration which is 
7- layer dielectric sheets. The wall is constructed by 3-sheets, 
namely 2- skins and one core. The radar target is high density 
dielectric or metal sheet separated 1.5m from the wall. Free 
space (air) is assigned to 1st, 5 l and 7 sheet. Figure 6 is 
reflection and transmission coefficients calculated by Equation 
(8) with no radar target (metal-like dielectric sheet) and Figure 
7 is with the sheet. The frequency range is 1-8 GHz, the 
incident angle is ~ 60 degree and the vertical axis is the 
calculated coefficient including its loss causing by complex 
permittivity (relative permeability jx r = 1). Difference between 
Figure 6 and 7 is due to the high density sheet, and frequency 
response in low band especially is changed and rippled signs 
appear according with incident angle. 

Generally, most of actual objects are made from various 
metallic and/or non-metallic materials. Consequently the 
general and systematic code is indeed necessary in order to 
predict characteristics of target RCS. According this novel idea 
of multiple reflections and these non-metallic multi-sheets, one 
can estimate scattering fields by complex objects with 
refracting wave like radome. It may be available for microwave 
remote-sensing application as fully polarimetric components 
regarding objects such as sea or ground as a single sheet with 
deep thickness. These coefficients for multiple layers can be 
also applied to design radome or radio absorption material. 

V. ADAPTIVE SUBDIVISION SCHEME 

The object surface of computational model consists of 
triangular facets. When shadowing or multiple reflections 
occurs on the triangular facets, only a part of the triangular 
facets that being illuminated by the incident or reflected wave 
have a contribution to the scattering field. We propose an 
adaptive subdivision new scheme to cope with this situation. 

At each specific incident angle or reflection angle, the 
coordinate of the triangular facets is transferred to a new 
coordinate that the z-coordinate is along with the direction of 
the incident or reflected electric wave. Then, all triangular 
facets of the model are projected onto the new xy-plane. The 
overlapped part between the current triangular facet and the 
projected facet is subtracted from the current facet and the 
current facet is subdivided into 3 triangular new facets in this 
case as shown in Figure 8. The contribution to the scattering 
field by the current triangular facet is re-organized as the 
contribution by the 3 newly generated triangular facets. Each 
engine, namely ISE, RSE, SSE and MRE, employs the 
adaptive subdivision scheme to cope with the shadowing. 

Shadowing example for a square trihedral coiner with a 
sphere is shown in Figure 9 and 10. The former is a case of no 
processing and the latter is shadow processing current intensity 
distribution at incident angle 8 =55,0 = 45 degrees. The 
contribution to the scattering field by the current triangular 
facet is reorganized as the contribution by the three triangular 
facets. This adaptive subdivision processing is applied to the 
each engine including Scattering Shadowing Engine and 
Multiple Reflection Engine of the developed code. 

In "Configure/Calculation" menu in the ACES opening 
window, the GUI consists of five modules: Mono-static RCS 
versus angle, Mono-static RCS versus frequency, Bi-static 
RCS versus angle, Bi-static RCS versus frequency and 
Induced current. In Figure 11, for reference, the simulation 
parameters of the developed code for mono-static RCS 
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calculation are displayed. The incident wave operates at 
frequency of 10GHz and its polarization is in ^-direction. For 
simple the material of the calculated model and the ground is 
PEC (perfecting electrical conductor). Incident Shadowing and 
Multiple Reflection are also considered. Figure 12 is outline of 
the developed code configuration and calculating flow. 



developed software code, which is available to calculate and 
estimate complex and large radar target with/out wide area 
like ground or sea. The code can also calculate fullpolarimetric 
fields. Using this code, RCS of corner reflector already verified 
to moment method and they coincide with each other within 
±0.5 dB [9]. 





Figure 8: Adaptive shadowing in subdivision processing. 
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Figure 9: Surf ace current distribution without shadowingon asquaretrihedral 
corner with a sphere (incident angle0 o = 55,0 = 45 deg). 
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Figure 10: Surface current distribution with shadowing on a square trihedral 
corner with a sphere (incident angle0 o = 55,0 = 45 deg). 

VI. SOME SIMULATION RESULTSTO RCS AND DOPPLER 
ESUMALITION 

This section presents a RCS simulation example of airplane 
as shown in Figure 13 and a large tree on wide area ground, 
and Doppler spectrum simulation of a helicopter using the 
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Figure 11: Mono-static RCS calculation parameters in the developed code. 
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Figure 12: Outline of ACES configuration and calculating flow. 

As shown in Figure 13, for the airplane, the x-, y-, and z- 
dimension are 9.9m, 6.0m, and 2.0m, respectively, and z- 
coordinate range is 10m. Figure 14 shows the polar plot of 
monostatic RCS at 10GHz without ground plane. We can see 
from the Figure that cross-polarized component of RCS is 
comparatively high level in a half space, 180 ~ 360 degree 
(rear area of the airplane) because the airplane has twisting 
facets to incident field vector direction. 

In radar remote-sensing field, analysis of SAR image by 
full polarimetric information is now in progress [10]. Image of 
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SAR or Inverse SAR can be regarded as spatial infinitesimal 
variation of the electromagnetic field, so SAR simulation and 
general radar signal processing like clutter suppression or 
Doppler spectrum become to be available by this code. 

Figure 15 and 16 show surface current distribution on a 
large tree (h = 10m) illuminated by 10GHz plane wave and 
mono-static co/cross-polar pattern which are 0(V)-polarized, 
respectively. The ground (50 x 50m 2 , s r = 5 + y'0.5,/i r = 1) 
with the electrically large tree (h » 3,333A) is non-metallic 
model. Figure 16 is co/cross-polarized RCS patterns for 
horizontal ground plane whose reflection coefficient is 
calculated by taking the thickness of first sheet over h » 
1,000/1 according to the previous discussion and we can find 
cross-polarized field components around +60 ~ +80 degree 
according to the tree branch. If off-nadir angle of SAR is 
around 30 degree, this model RCS, a, is a value at 8 = +30 
degree. Additionally speaking, there is need in future, to 
simulate taking account of normalized RCS per unit area or 
beam width of SAR antenna, surface terrain, shape and 
permittivity of trees, and some verification with other code or 
experimental comparison etc. 
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Figure 13: Airplane model developed in the ACES Code. 
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Figure 14: Monostatic RCS polar plot at 10GHz, polarized W (09) and cross- 
polarized VH (9(p) for the airplane model in Fig. 12. 
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Figure 15: Tree model on a wide ground and surface current. 




Figure 16: Co/cross-polarized RCS pattern forthe model in Figure 15. 

Lastly we show a dynamic simulation example for Doppler 
estimation. Figure 17 is a model for helicopter with 4-main 
rotors to use in the developed code. The diameter of main and 
tail rotor is 16.4m and 4.7m, respectively, and rotational speed 
at its edge is 216 and 121 m/s. Then each rotational number 
and periodic cycle is 4.2 round/s and 8.4 round/s, and 14.4 kHz 
and 8.1 kHz. If PRF (pulse repetition frequency) of radar 
transmitted pulse is 50 kHz, sampling of this model simulation 
should be fitted to this PRF. In this situation, rotating main and 
tail rotors and calculating of the whole helicopter under 
condition of hovering in the code, we obtain simulated results 
shown in Figure 18 to 20. Figure 18 is RCS change with 
respect to time (rotation). The RCS amplitude around center is 
not so changed and caused by body area of the helicopter and 
high level amplitude is by 4-main rotors which peak are 
occurred per 0.24/4s because of 4-rotors 1/4.2=0. 24s. Figure 19 
is spectrum processed by short time FFT (Fast Fourier 
Transformation) of data in Figure 18. Figure 20 also is 
processed in short time (512 points). From Figure 20, we can 
see locus of the main rotor with maximum level 14.4 kHz 
which is mentioned as the round cycle of the main rotor. 

Furthermore in order to verify the above simulation we 
measured Doppler spectrum using X-band pulse Doppler radar. 
The radar target is an actual helicopter which shape is like the 
model in Figure 17. Figure 21 and 22 are processed data of the 
measurement as same parameters without incident angle in 
Figure 19 and 20, respectively. This elevation angle from a 
ground at the measurement is 20 ~ 30 degree. This difference 
of incident angle affect to the mismatch around zero-Doppler 
spectrum of these Figures. And, unfortunately, non-linear 
devices like a RF amplifier in the radar are not taken into 
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account to the simulation and IF band of the radar receiver is 
within 15 kHz because of its low pass filter. 
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Figure 17: Helicopter simulation model with 4main rotors. 



E 

■-T, 









MJt 


i- 


■•••!•• 


: 


..:. ..!..... 




.. .. 


14 .ft 










! 


'. 


'. 





(I.l 



(1.2 0.3 

Time(s) 



Figure 18: Time varianceof simulated RCS at 10 GHz at normal incidence 
to nose of the model in Figure 17. 
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Figure 19: SimulatedDoppler spectrum using Figure 18. 
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Figure 20: Simulated Short -time spectrum usingFigure 18. 
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Figure 21: Measured Doppler spectrum by using the pulse Doppler radar: 
Range /div is different from Figure 19. 
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Figure 22: Measured short -time spectrum usingFigure 21. 



VII. CONCLUSION 

In this paper, consciously considering large and complex 
object with non-metal portion and semi-infinite space like sea 
or ground, we proposed the systematic RCS prediction code 
using accurate high-frequency asymptotic methods, PO and 
PTD which are based on electromagnetic surface and edge 
currents. Moreover the code can calculate cross-polarized field 
for full polarimetric analysis of radar remote-sensing imaging 
and decomposition. For complex shape objects, software code 
has to treat not only first reflection but also secondary 
reflection and diffraction. From this reason, algorithm of 
shadowing is always necessary. For full polarized RCS, the 
developed code has an optimum processing unit engine 
adaptively to treat re-meshing, recurrently calculating unit for 
multi-layered dielectric sheets and multiple scattering which 
described in this paper. The trial version of this developed code 
is available to download from URL [11]. 

Additionally, we indicated a dynamic simulation example, 
Doppler analysis using this developed code, and then verified it 
with measurement data of an actual helicopter by the Doppler 
radar. These Doppler processing and analysis technology can 
directly be applied to NCTR (Non Co-operative Target 
Recognition). 
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